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Abstract 



The multiphoton ionization and photoelectron spectra of magnesium were 
studied at laser intensities of up to 6 x 10^'^ Wcm~^ using 150 fs laser pulses 
of a wavelength of 400 nm. The results indicated that a variety of different 
ionization mechanisms played a role in both types of spectra. A theoretical 
model describing the processes is presented and the routes to ionization are 
identified. The work demonstrates the significance of the two-electron nature 
of the atom in interpreting the experimental results. 
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I. INTRODUCTION 



The behavior of atoms in intense electro-magnetic fields has been under investigation 
since the advent of the laser. Especially since the development of short-pulse, high-intensity 
systems, experiments have revealed a multitude of novel effects such as multiphoton ioniza- 
tion (MPI) above threshold ionization (ATI) p|, harmonic generation etc. Most of 
the theoretical models originally proposed to describe the observations involved essentially 
a single electron interacting with the electromagnetic field. This approximation has led to 
significant understanding of the atomic behavior in the field. However, the multi-electron 
nature of most atoms can play a significant role under certain circumstances. Theory then 
has to take into account multi-electron excitation, correlations between the electrons, mul- 
tiple continua, i.e. the detailed atomic structure in order to be able to predict the atomic 
behavior. This fact emerges as both theoretical and experimental techniques become more 
advanced, enabling the exploration of new regimes and the detailed testing of theoretical 
models. 

In multi-electron atoms, significant double excitation can take place leaving the ion not 
in the ionic ground state but in an excited state or possibly ejecting a second electron. 
In the alkaline earth atoms for example, the existence of doubly excited-state manifolds 
between the two ionic thresholds alters drastically the single-electron picture rendering it 
unsuitable. The conditions under which these effects may occur are easily within reach 
of present experimental setups. Past work by Dalwoo Kim et al [^], DiMauro et al 
and van Druten et al |^ has looked into the multiphoton ionization of magnesium with 
nanosecond and picosecond laser pulse durations and indicated the importance of short 
pulse experiments in understanding the processes involved. At the same time, theoretical 
methods incorporating these effects are also being developed creating the need for 

further experimental investigations. 

In the present work we seek to investigate the extent to which electron correlations play a 
role in the high intensity regime for multi-electron atoms. We chose magnesium, a fairly low 
atomic weight alkaline earth element, as it is accessible to both experimental and theoretical 
investigation. We report measurements of the multiphoton single and double ionization and 
ATI of magnesium using 150 fs, 400 nm laser pulses. A theoretical model has been developed 
and applied which describes the processes leading to ionization and shows that the inclusion 
of one and two-electron effects are necessary for an accurate description. Measurements of 
the multiphoton ionization yield displayed the expected perturbative behavior for the pro- 
duction of the Mg"*" ion but showed unexpected features in the double ionization yield. The 
ATI electron spectra also exhibited complex structure, originating from a variety of single 
and double ionization mechanisms, whose relative importance changes as the intensity is 
varied. The theoretical treatment involved both a perturbation theory description including 
all correlation effects between the active electrons and additionally, for the ionization of 
Mg"*", a non-perturbative approach. 



II. METHOD 
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A. Experimental set-up 



A plan of the experimental set-up is shown in figure |I|. The laser employed consists 
of a titanium sapphire oscillator amplified in a regenerative amplifier utilizing a chirped 
pulse amplification technique [0. The laser operates at 1 kHz repetition rate, with a 
measured pulse duration of :^150 fs at a wavelength of ~798 nm. The second harmonic of 
this beam was then produced using a focused geometry in a 250 fim thick BBO crystal. The 
harmonic radiation was separated from the co-propagating fundamental by employing three 
dichroic mirrors that were transparent to the fundamental (denoted by M2i^ in figure 0). The 
harmonic generation efficiency of this setup was ~10%. The laser energy of the fundamental 
was set to be less than 300 fiJ, thus producing a maximum of ~30 yuJ of second harmonic 
radiation. 

The resulting harmonic laser beam was focused in a time-of-fiight spectrometer (TOF) 
by a MgF lens of 15 cm focal length. A half-wave plate was positioned prior to that in 
order to align the polarization of the light with the TOF tube axis. Under the conditions 
of the experiment, the lens did not affect significantly the beam either in terms of non- 



linear absorption or by pulse- front curvature at the focus |irT|. The energy measurement 



was achieved with a fast vacuum photodiode which was cross-calibrated with an average 
power meter. 

The time-of-fiight spectrometer, which has a field-free fiight tube of 15 cm oriented 
perpendicular to the laser propagation axis, is equipped with a //-metal magnetic shield. 
This shield, in combination with a repeller plate and the appropriate polarity of the charged 
particle detector, allowed the use of the instrument alternately as an electron spectrometer 
and as an ion time-of-fiight system. For the electron energy measurements the repeller was 
grounded, whereas for the operation of the TOF as a mass spectrometer the repeller was 
charged with a voltage of up to +1000 V. The charged particle detector at the end of the 
fiight tube was a double micro-channel plate in a chevron configuration. Its front side was 
grounded for detecting electrons and at a negative voltage for ions. 

An oven was inserted at the top of the interaction chamber filled with magnesium pellets 
which, when heated, produced an atomic beam as the vapor exited an aperture. The atomic 
beam met the laser beam and the time-of-fiight tube at right angles. 

The acquisition system was essentially the same for both electron and ion detection. The 
signal from the micro-channel plates first went through a fast amplifier and a discriminator. 
Subsequently, its time of arrival with respect to the laser pulse was recorded by a multiple- 
event time-digitizer card which was installed in a personal computer. The main feature of 
this card was that it could acquire and time one event per 0.5 ns time-bin with no dead- 
time between bins. Combined with the bandwidth of the cabling and electronic processing 
stages, the effective time resolution of the system was approximately one nanosecond. This 
acquisition system was not sensitive to the amplitude of the signal and hence it was operated 
in the event-counting regime. This meant that the number of events in each bin had to be 
less than one per laser pulse, a condition which was ensured throughout the measurement. 

At the higher laser intensities, for which the counting method was no longer appropriate 
due to the increased ionization rates, the analogue signal was obtained directly from the 
micro-channel plates and recorded on a digital oscilloscope. For the ion yield versus laser 
intensity spectra, where data were recorded using both methods of acquisition, the data 
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were matched at several points in an intermediate regime. They were further cross-checked 
with data from the less abundant isotopes of magnesium which would be still in the counting 
regime at intensities for which the main peak was not. Lastly, to ensure that no space-charge 
effects were distorting the data, measurements were taken at different oven temperatures 
(and consequently atomic number densities) and compared. No rigorous absolute energy 
calibration of the electron spectrometer was performed. The calibration of the energy scale 
was based on the expected inter-peak separation of the several ATI peaks. This resulted 
in a determination of the relative peak positions to better than 0.1 eV, although the de- 
termination of the absolute peak positions was hindered by uncertainty in the zero time. 
Comparison to data with xenon was consistent with the magnesium results in so far as the 
energy calibration of the spectra was concerned, confirming the validity of the calibration. 
Xenon was also used to estimate the intensity of the laser at the interaction region [p!^ -P 
and compare it with the number obtained by estimating the focal spot size. 



B. Theoretical modeling 



We now propose a theory to describe the experimental results and explain most of the 
essential features. Magnesium in the present context is a rather complicated system since 
it involves many coupled states and different ionization channels. However, in that case, 
we were able to identify and isolate atomic states and transitions connected to the leading 
processes. These elements are summarized in figure ^. We consider that the electrons mea- 
sured in the experiment are produced by 4 different processes. Two lead to the production 
of Mg"*" (single ionization to the ground state Mg"*"(3s) by absorption of 3 photons or to the 
first exited state Mg+(3p) by absorption of 4 photons) and are labeled process (1) and (2) 
respectively. Two other processes sequentially double ionize the system from the ion Mg"*". 
Mg"*""*" is thus produced by absorption of 4 photons from Mg"'"(3p) or by absorption of 5 
photons from Mg"'"(3s). The basic assumption in the theoretical approach is that double 
ionization is produced through a sequential process. Therefore, double ionization through 
direct (coherent) absorption of 8 photons (see path (5) in figure 0) is neglected. Under these 
conditions, the population of each species is obtained by solving the following rate equations: 

^3. = r£L 3s., N,,. - rtl^,^,^ Ns. (1) 

with the initial conditions N^s'^lt = — oo) = 1 and all other populations set to zero. A'3^2, 
and N^p represent the time-dependent populations of Mg(3s^), Mg"'"(3s) and Mg"'"(3p) 
respectively. The time-dependent width Tg-^^^^^^^^ refers to the p-photon transition from the 
initial state 3/ni to the final state ^271^3. When the specification of the angular momentum has 
been removed from the state notation (like in ^2^3), the width includes the contribution from 
each accessible symmetry, kn refers to the momentum of the electron ejected through the 
process labeled (n) in figure ^ The electron energies (reported in table ||) are deduced from 
the energy conservation principle excluding laser induced level shifts. The atomic structure 
calculations providing the atomic energy levels and the transition cross-sections have been 
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described in detail in a precedent paper . The energies and wave-functions are computed 
using a frozen-core configuration interaction (CI) procedure. The core Mg"'""'"(ls^2s^2p^) 
(where electrons are assumed "frozen" ) is represented by a self-consistent-field wave function. 
The core potential (including a dielectronic polarization term) is first used to determine 
Mg"*" orbitals. We then compute two-electron wave-functions by diagonalizing the 2-electron 
Hamiltonian in a basis of Mg"*" orbitals configurations. The multiphoton ionization rates 
r3s2^3;fc ^'^d T^^il^k'k ^-^s calculated within lowest order perturbation theory (for the 



calculation of MPI matrix elements see also [|T7| and references therein). The rates leading 



to double ionization from Mg"*" are deduced from non-perturbative ATI calculations [T^ (see 
below) and are in agreement with the perturbative results. The spatio-temporal dependence 
of the laser is accounted in this model by assuming the following field envelope: 

where we assume a cylindrical symmetry (r, z) for the space macro dependence of the in- 
tensity. K and r are respectively the spatial and temporal full width at half maximum 
(FWHM). 

Besides the perturbative evaluations of ionization rates, we have performed non- 
perturbative calculations of ionization from Mg"'"(3s) and Mg"'"(3p) by numerically solving 
the time dependent Schrodinger equation (TDSE) [|^. A pseudo-potential JlQl is used to 



model the Mg"*" ionic potential. The total wave-function is expanded in terms of spherical 
harmonics and radial B-spline functions. The total Hamiltonian is then propagated in time 
and the electron spectrum is calculated at the end of the pulse. This approach has the 
advantage of neglecting none of the couplings. It therefore accounts for all dynamical level 
shifts and possibly resonant transitions. These latter effects are of importance to interpret 
some of the experimental results as it is discussed below. 



III. EXPERIMENTAL RESULTS AND COMPARISON WITH THEORY 

A. Ion Yields 

Figure ^ shows a typical time-of-fiight ion mass spectrum obtained in the experiment, 
in logarithmic scale. The two charge states of magnesium (Mg+ and Mg"'"^) are evident, 
each displaying peaks corresponding to its three dominant isotopes. The isotope ratios are 
as expected and no presence of impurities is detected. By integrating under the main peaks 
of such spectra obtained for different laser energies, the ion yield curves of singly and doubly 
ionized magnesium were obtained. 

The resulting ion yield curves are shown in figure ^. The theoretical curves (given by 
solving equation |l| and denoted by the solid lines) agree very well with experimental values 
over most of the range investigated. Small discrepancies in the heavily saturated region 
above 2 x 10^^ W cm~^ are expected since the detection volume was different to the inter- 
action volume. The precise experimental intensity was determined by the agreement with 
the theoretical curves. However, some uncertainty to this estimate can be expected due to 
the fact that the macro- dependence of the field envelope assumed in the theory strongly 
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affects the ion yield behavior, especially in the saturated region, sometimes even affecting 
the apparent saturation intensity. The theoretical rate for three photon ionization of Mg 
is r^^l^3^^^ = 0.45 X 10*^ (j^)^ a.u. (where Iq = 6.44 x 10^^ W cm"^). The corresponding 

generalized cross-section o'^j2^^gi^^ = 0.85 x 10~^° cm^ sec^ is in agreement with the values 
calculated by Chang and Tang One can easily calculate that single ionization satu- 

rates at about / = 5 x 10^^ W cm^ for a pulse duration of 150 fs, in agreement with the 
experimental ion yield shown in figure |^. 

In the lower energy range of the figure (around I = 4x 10^^ W cm~^) a clear enhancement 
to the Mg^^ yield relative to what expected from the theoretical results can be seen. This 
enhancement seems to level off, the experimental curve tending back to the theoretical one, 
at the onset of saturation of single ionization. This may be an indication that two-electron 
ejection is connected to the population of the atomic ground state which would be the case 
for e.g. non-sequential double ionization. It should be noted that the Keldysh or adiabaticity 
parameter 7 = ^Ip/{2Up) (where Ip is the ionization potential of the atom and Up ~ 
9.33 X 10~^^/(W cm~^)A^(/im) is the ponderomotive energy of a free electron in a laser field 
of wavelength A and intensity J) is not favorable to observe tunneling, rescattering or shake- 
off effects (7 = 2 at the highest intensity). It is therefore possible that the enhancement is 
related to atomic structure, or to dynamical shifts leading to either resonances or channel 
closing as it is discussed in section |111 B 1| . 

B. Electron results 

A typical electron energy spectrum is shown in figure |^. The spectrum displays a set 
of main ATI peaks, each accompanied by a pronounced side peak (at approximately 1.3 eV 
to its low energy side). While the main ATI peak can safely be attributed to three-photon 
ionization of neutral magnesium leading to the ion Mg"'"(3s), the origin of the side peak is 
not so clear without further investigation. Note that the main peaks in figure ^ are offset 
compared to the expected energies corresponding to normal perturbation theory. This is an 
experimental artifact due to the lack of absolute energy calibration of the TOF as it is not 
relevant to the present investigation (see also section p.1 A] ). 

Table H gives the calculated energy of the ejected electrons upon ionization through 
the different processes depicted in figure ^. Comparing these with the electron spectra of 
figure ^ suggests two possible mechanisms that would result in the side peak in the recorded 
data: (i) 4-photon ionization of Mg leading to the Mg"'"(3p) threshold, i.e. the first excited 
state of the magnesium ion (process (2) in table |). That would give an energy difference 
between the main and side peak of (1.656 — 0.34) = 1.316 eV, (ii) 5-photon ionization of 
Mg"'"(3s) leading to the 2p^ ground state of Mg++ ((4) in table |), which would give an 
energy difference of (1.656 — 0.47) = 1.186 eV. Both of these candidates are within the 
experimental uncertainty and the electron peak widths, especially at the higher intensities. 
Note that direct double ionization (process (5) in table |I|) would give a continuous electron 
distribution with maximum energy of 2.13 eV, as that energy is shared between the two 
electrons. 
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1. Electron shifts 



Careful analysis of the electron spectra can give more insight into the mechanisms in- 
volved. Figure || shows a detail of the ATI spectra of figure | with additional plots for 
different intensities. If one concentrates just on the peak positions one notices that the main 
peak shifts smoothly up to a point and then stops completely. In contrast, the side peak 
initially shifts, then stops and finally starts shifting again. This behavior is shown clearly 
in figure |^ where the positions of two electron peaks are plotted versus laser intensity. The 
square points correspond to the main peak while the circles to the side peak. 

The main peak originates from the processes (1) and (3) shown in figure ^ Our calcu- 
lations show that, in all cases, the contribution from process (1) (three-photon ionization 
leading to a Mg~''(3s) ion) dominates over (3). Indeed, this latter process originates from 
Mg"'"(3p), whose production from Mg(3s^) requires the absorption of 4 photons in contrast 
to the 3 photons involved in process (1). The flattening can be explained in terms of the 



ion production saturating p!^J13|| . As the intensity increases, there comes a point where, in 
effect, the atom will ionize before the end of the laser pulse, and so it will not experience the 
full power of the laser. Therefore the curve levels out. To account for the specific, detailed 
behavior one has to take into account all the factors that would contribute to the shift, e.g. 
(i) ponderomotive shifts and (ii) AC Stark shifts. 

The lower data points in figure (corresponding to the side peak) have a different be- 
havior, which cannot be explained solely by saturation. The two different trends observable 
in the energy shift indicate that more than one process is at work. As can be seen from the 
theoretical calculations in figure |Ka) and (b), the side peak is a (unresolved in our data) 
combination of two processes, one coming from single ionization (2) the other from double 
sequential ionization (4) (see figure^). At low intensity (figure ||(a)), the calculations show 
that the side peak is dominated by single ionization (2). At about / = 1.8 x 10^^ W cm^^, 
processes (2) and (4) have an equal contribution, but for higher intensities the process (4) 
largely dominates over (2) (figure |^(b)). Figure |](c) and (d) indicates the development of 
the relevant atomic populations as the laser pulse evolves. In fact, as we will see below, 
the process (2) should have a significantly smaller contribution for intensities larger than 
/ = 8 X 10^^ W cm~^. The electron production originating from the decay to the first excited 
ionic state would saturate at the same intensity as the electrons of the main peak, in agree- 
ment with the first trend observed in the data. On the other hand, the electrons originating 
from sequential double ionization (4) are expected to saturate at higher intensities. As the 
intensity increases, this second process dominates in the side peak, and the change of slope 
observed after the initial leveling off is due to the shift experienced by this second different 
peak. This latter shift is shown in figure |^ where we plot the energy shift of electrons orig- 
inating from process (4) versus the intensity. This graph results from a non-perturbative 
calculation which includes both AC Stark shift and ponderomotive shift. Note that the shift 
is not uniquely due to the ponderomotive effect as is often the case in ATI measurements. 
This is shown in figure ^ where we detail the contribution of the ponderomotive shift to the 
total shift. It should be noted that the experimental shifts appear to have a smaller magni- 
tude than the theoretical ones, indicating a smaller experimental intensity than previously 
extracted from the ion yields. However, considering that the acquisition of the electron 
signal takes place under different conditions to that of the ion signal (field-free, as opposed 
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to a high apphed voltage for the ion extraction), it is extremely likely that the extraction 
region, and therefore the effective intensity seen by each species could well be different. The 
uncertainty in the determination of the intensity in the ion yields should also be taken into 
account, as described in a previous section. 

Our calculations show that, for intensities larger than 10^^ W cm~^, the double ionization 
is dominated by the process (1+4) (see figure At lower intensities, the channel (2+3) 
has a significant contribution to the production of Mg++. As already mentioned, the double 
ionization yield displays a significant deviation at low intensities. We considered the laser 
induced energy shift of the Rydberg series Mg{3pnl) (converging, as — > oo, to the sum 
of the AC Stark shift of Mg+(3p) and the ponderomotive shift Up) and the AC Stark shift 
of Mg(3s^). Through this we found that, for intensities larger than 8 x 10^^ W cm~^, 
the 4-photon ionization channel (2) is closed and ionization to Mg"'"(3p) now requires the 
absorption of 5 photons. The effect is that the experimental double ionization yield deviates 
from the normal picture at about / = 8 x 10^^ W cm~^. Note that the single ionization 
channel, dominated by the process (1), is almost not affected by this effect. Nevertheless, 
we see in figure ^ that the theoretical Mg"'"^ yield lies below the experimental curve in 
the low intensity region. According to the previous considerations, it should fit the Mg^^ 
experimental yield at low intensity and overestimate it for / > 8 x 10^^ W cm~^ (in the rate 
equation system (|^), 4-photon ionization of Mg to Mg+(3p) is allowed at all intensities). 
Thus the discrepancy between theoretical and experimental results in the latter region still 
remains to be explained. Other effects, like direct and/or sequential ionization (eventually 
involving higher thresholds) enhanced by the presence of resonances, should be explored. 
Consequently, the sequential model adopted here for double ionization may not be valid at 
low intensities, where further experimental and theoretical efforts may be required in order 
to complete the present analysis. 

IV. CONCLUSIONS 

The results presented in this paper are some of the very few works studying in both 
experimental and theoretical detail the MPI and ATI spectra of an alkaline earth using ultra- 
short pulses. The work demonstrates clearly the significance of the two-electron nature of 
the atom in these spectra, both in terms of the role of the electron correlations, and with 
respect to the precise level structure of both the atom and the ion. 

The attribution of the pronounced sets of side peaks, accompanying the usual ATI spec- 
trum, was not trivial, and required a theoretical investigation taking into account all of the 
above factors in order to explain their origin. As it turned out, more than one processes 
contributed to their appearance, each dominating at a different range of intensities. Most 
notably, one of the two dominant processes at low intensity was the decay to the first excited 
electronic state of the ion, a pathway which involves the excitation of both valence electrons 
and coupled ionization channels. However, the prominence of this process was diminished 
as the intensity increased due to the depletion of the ground atomic state, and overtaken 
by the electron originating from sequential double ionization from the ground state of the 
atom. In principle, another new peak should arise due to sequential double ionization from 
the excited ionic state, but the position of this was calculated to be too close to the main 
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structure (thrcc-photon ionization of Mg) to allow for its observation. It might be useful to 
conduct an extension of this stTidy to a different photon energy, at the moment not available 
in our laboratory, at which these proximities are avoided, and at which the peaks from each 
of the processes are well resolved. 

The detailed studying of the electron energy shifts of both the main and the side set 
of peaks also pointed to the importance of the factors mentioned above. Quite apart from 
their being used as a tool towards the identification of the side peak, the shifts displayed a 
behavior not sufficiently explained by the common analysis in terms of ponderomotive shifts. 
Due to the complex level structure of both the atom and the ion, with transitions close to 
resonance with the laser field, AC Stark shifts were found to play an important role and could 
no longer be ignored. The multiphoton ionization yields, despite the relatively low intensities 
employed, also display an interesting enhancement of double ionization which deviates from 
the normal picture. We have found that, due to AC and ponderomotive shifts, the channel 
(2) requires the absorption of 5 photons for intensities larger than 8 x lO^^W cm~^. This 
is a possible explanation for the enhancement of double photon ionization, but other effects 
should be explored, in particular the possibility of near-resonant enhancement or, more 
tentatively, direct double ionization, possibly enhanced by electron correlations and presence 
of multiple doubly excited manifolds. 

All of the above indicate that a full picture of the atom needs to be explicitly contained in 
any investigation attempting to elucidate the high-intensity behavior of two electron atoms, 
and possibly other atoms as well. 
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FIG. 1. Experimental setup. SHG: Second harmonic generation, M.2uj: mirror for the second 
harmonic radiation, TOF: Time-of-flight, MCP: Micro-channel plates. 
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FIG. 2. Partial energy level diagram of magnesium and the relevant transitions involving 

400 nm photons in the neutral atom (Mg) and its ionic states (Mg+ and Mg'^^). The different 
paths to ionization indicated by numbers (1) to (5) are described in the main text and in Table |. 
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FIG. 3. Typical time of flight spectrum obtained during the experiment. While the singly and 
doubly ionized magnesium species can be observed, no impurities are evident. The isotopic ratios 
are as expected (A=24: 78%, 25: 10%, 26: 11%). Note the logarithmic scale. 
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FIG. 4. Experimental and theoretical ion yields. The square data points correspond to 
Mg"*" whereas the circles to Mg"*""*". The solid lines are the calculated yields (rate equations and 
perturbation theory). 




2 4 6 8 10 12 14 

Energy (eV) 



FIG. 5. Electron spectrum at a laser intensity of 2 x lO^^Wcm ^. Two dominant sets of 
electron peaks are observable. 
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FIG. 6. Enlarged section of the electron energy spectra of figure |5| focusing on a main peak 
and a sub-peak, for various laser intensities. A shift of the peaks towards the higher energies as 
the laser power decreases is observable. (Intensities from top to bottom in 10^^ W cm~^: 4.7, 3.6, 
2.3, 1.9, 1.4, 0.7, 0.4, 0.24, 0.22). 
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FIG. 7. Electron peak positions vs. laser intensity for a pair of peaks at approximately 8.5 eV 
and 7.1 eV. Note the axis break, (see text for more details). 




FIG. 8. Ion and electron populations calculated with no intensity space macro-dependence at 
two different intensities, (a), (b) Electron production at intensities of / = 1.7 x 10^^ W cm~^ and 
/ = 3 X 10^^ W cm~^ respectively. Electron peaks correspond to processes 1-4, as indicated in figure 
(c), (d) Time dependent populations of different ionic species at the intensities corresponding to 
cases (a) and (b) as the laser pulse evolves, (solid line): Laser intensity normalized to the maximum 
intensity; (dashed line): ground state of the atom Mg(3s^) (short-dash line): ground state of the 
ion Mg"'"(3s) produced by three-photon ionization from Mg(3s^) along with electron of energy El; 
(dotted line): first excited state of the ion Mg"'"(3p) produced by four-photon ionization from 
Mg(3s^) along with electron of energy E2; (dash-dot-dot line): doubly ionized Mg"'""'"(2p^) through 
four-photon ionization of Mg"'"(3p) along with electron of energy E3; (dash-dot line): doubly ionized 
Mg^+(2p^) through five-photon ionization of Mg^(3s) along with electron of energy E4; 
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FIG. 9. Total shift (solid line) in energy of electrons originating from process (4), resulting 
from a non-perturbative calculation of the ionization of Mg+(3s). Also shown (dashed line) is the 
contribution to the total due to the Mg+(3s) AC stark shift. 



TABLE I. Table of the energies of the different ionic thresholds in which the ionization may 
leave the Mg atoms, the number of 3.1 eV photons that need to be absorbed to reach each of 
these (-/Vph) and the energy of the ejected electron upon ionization (£^ei)(All numbers refer to field 
free energies, with energy shifts not taken into account). The various processes are illustrated in 
figure |2|. 



No Transition Energy (eV) A^ph £'ei(eV) 

(T) Mg(3s2)-Mg+(3s) 7.644 3 1.656 

(2) Mg(3s2)-Mg+(3p) 12.06 4 0.34 

(3) Mg+(3p)-Mg++(2p6) 10.614 4 1.786 

(4) Mg+(3s)-Mg++(2p6) 15.03 5 0.47 

(5) Mg(3s2)-Mg++(2p6) 22.674 8 2.13 
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